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ABSTRACT 

The Wide-field Infrared Survey Explorer (WISE) has completed its all-sky 
survey in four channels at 3.4-22 /im, detecting hundreds of millions of galaxies. 
We merge the WISE mid-infrared data with optical data from the Sloan Digital 
Sky Survey (SDSS) and provide a phenomenological characterization of WISE 
extragalactic sources. WISE is most sensitive at 3.4 /xm (Wl) and least sensitive 
at 22 /im (W4). The Wl band probes massive early-type galaxies out to z > 1. 
This is more distant than SDSS identifies early-type galaxies, consistent with the 
fact that 28% of 3.4yum sources have faint or no r-band counterparts (r > 22.2). 
In contrast, 92 - 95% of 12 /im and 22 fim sources have SDSS optical counterparts 
with r < 22.2. WISE 3.4 /im detects 89.8% of the entire SDSS QSO catalog at 
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SNRiy! >7a, but only 18.9% at 22 fjm with SNR W4 >5a. We show that WISE 
colors alone are effective in isolating stars (or local early- type galaxies), star- 
forming galaxies and strong AGN/QSOs at z < 3. We highlight three major 
applications of WISE colors: (1) Selection of strong AGN/QSOs at z < 3 using 
Wl — W2 > 0.8 and W2 < 15.2 criteria, producing a better census of this 
population. The surface density of these strong AGN/QSO candidates is 67.5 ± 
0.14 per deg 2 . (2) Selection of dust-obscured, type-2 AGN/QSO candidates. We 
show that WISE Wl-W2> 0.8, W2 < 15.2 combined with r - W2 > 6 (Vega) 
colors can be used to identify type-2 AGN candidates. The fraction of these type- 
2 AGN candidates is l/3rd of all WISE color-selected AGNs. (3) Selection of 
ultra-luminous infrared galaxies at z ~ 2 with extremely red colors, r — WA > 14 
or well-detected 22 /zm sources lacking detections in the 3.4 and 4.6 jum bands. 
The surface density of z ~ 2 ULIRG candidates selected with r — W4 > 14 
is 0.9 ± 0.07 per deg 2 at SNRvk4 > 5 (the corresponding, lowest flux density of 
2.5 mJy), which is consistent with that inferred from smaller area Spitzer surveys. 
Optical spectroscopy of a small number of these high-redshift ULIRG candidates 
confirms our selection, and reveals a possible trend that optically fainter or r—WA 
redder candidates are at higher redshifts. 

Subject headings: galaxies: infrared luminous - galaxies: starburst - galaxies: 
high-redshift - galaxies: evolution 



1. Introduction 

With the power of large number statistics, wide-area astronomical surveys in this decade 
have made many significant discoveries. A new generation of high-precision, wide-area op- 
tical imaging surveys are being proposed for the end of this decade to address fundamental 
questions in cosmology, galaxy formation and to search for extrasolar planets. Most re- 
cently, the advent of the Wide-field Infrared Survey Explorer ( WISE] Wright et al.||2010" ) has 



provided the community with an unprecedented dataset in the mid-infrared. This unique 
NASA mission mapped the entire sky in four bands at 3.4, 4.6, 12 and 22 fim (Wl thru 
Wi), with 5(7 point source sensitivities better than 0.05, 0.1, 0.75 and 6 mJy, respectively. 
WISE 12 fim images are more than 100 times deeper than previous all-sky infrared survey 



missions, such as that provided by the Infrared Astronomical Satellite (IRAS; Neugebauer 



et~al]p84| ), while the 3.4 jj,m data is 1.5 magnitude deeper than the Two-Micron All-Sky 



Survey (2MASS; |Skrutskie et aLp006) K s data for sources with spectral energy distributions 



similar to an AO star; WISE is even more sensitive for red sources like K-stars, early-type 
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galaxies, AGN, and dust-obscured galaxies. 



WISE has made two public data releases, including processed atlas images and source 
catalogs. The first one was the preliminary release in April 2011, covering roughly half of 
the sky, and the second was the all-sky data release in March 2012, covering the entire sky. 
Although WISE has a significantly smaller aperture and the angular resolution is only half 
that of Spitzer, its unique all-sky coverage enables selection of large samples of extragalactic 
sources, allowing statistical studies of stellar photospheric emission at 3.4 and 4.6 /im and 
dust emission at 12 and 22 /im. The primary goal of this paper is to provide an empirical 
characterization of WISE extragalactic sources, and to identify the types of sources which 
can be isolated using mid-infrared colors. To achieve this goal, we combine WISE mid- 
infrared data with optical data from the Sloan Digital Sky Survey (SDSS) seventh data 



release (DR7; Abazajian et al. 2009). We characterize the color distributions and source 
types using broad-band photometry and spectroscopic information. The motivation for such 
a phenomenological study is to give the community a summary of the observational properties 
and the limitations of two large surveys, particularly for mid-infrared extragalactic objects. 
More detailed, quantitative analyses of specific extragalactic populations are discussed in 
companion papers by the WISE extragalactic science team. Specifically, Griffith et al. 



(2011) and Tsai et al. (2012a) discuss WISE selection of low-metallicity blue compact 



dwarf galaxies at z < 0.1. Donoso et al. (2012) examines the origin of 12 /xm emission in mid- 
infrared galaxies at z ~ 0.1 using WISE data in conjunction with SDSS DR7 spectroscopic 



data, while Lake et al. (2012) uses Keck spectroscopy to study the redshift distribution 



of flux-limited WISE samples. Stern et al. (2012) and Assef et al. (2012) present detailed 



studies of WISE-selected AGN within the COSMOS and Bootes fields, respectively, including 
careful analyses of completeness and reliability of WISE AGN selection using W1 — W2 color. 
Eisenhardt et al. (2012) presents the first results for z ~ 2 ultraluminous infrared galaxies 



(ULIRGs) discovered by WISE. Wu et al. (2012) and Yan et al. (2012) present the far- 



infrared properties and spectral energy distributions (SEDs) of similarly selected galaxies 
based on ground-based millimeter and Herschel space-based far-infrared data, respectively. 



Bridge et al. (2012) discusses the interesting Lya properties of the WISE ULIRG population, 



emphasizing the high rate of extended emission, so-called "Lya blobs". Tsai et al. (2012b) 
and Stern et al. (in prep) present detailed studies of two interesting H / 7S'i?-selected AGN, 



while Jarrett et al. ( |2012 ) and Petty et al. (2012) study spatially resolved, local (z < 0.1 



galaxies with WISE. Finally, Blain et al. (2012) presents the WISE properties of z > 6 



optically selected QSOs. And Gettings et al. (2012) discusses the first results of using WISE 



to identify high-redshift galaxy clusters. 

The organization of this paper is as follows. Section [2] describes the WISE and SDSS 
data, including the WISE sample selection and the SDSS optical data used by this paper. 
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Section [3] presents the main results, providing the technical details of the catalog matching 
and summarizing the brightness, color and photometric redshift distributions for galaxies 
detected by both WISE and SDSS (§ 3.1). We identify the characteristic color criteria which 
can be used to select large samples of strong AGN/QSOs (§ 3.2), type-2 AGN candidates 
(§ 3.3) and potential z ~ 2 ULIRGs (§ 3.4). Section [4] summarizes the main conclusions and 
discusses the implications for future studies using WISE. Throughout the paper, we adopt 
an Q M = 0.27, Q\ = 0.73 and H = 71 km s _1 Mpc -1 cosmology. 

2. Data 

2.1. WISE-SDSS sky coverage 

The data used by this paper are drawn from the overlapping sky region where both 
WISE preliminary public release data and SDSS DR7 data are available. Although this 
area is only l/4th of the total SDSS area, the number of galaxies included in our analysis 
is > 10 7 , large enough for the purposes of this paper. Because this paper started with the 
preliminary release data, we kept the same overlap region between the preliminary release 
data and SDSS DR7, but the actual WISE photometry are taken from the all-sky data 
release. Figure [T] presents the full-sky coverage map from WISE, with colors indicating the 
number of repeated single exposures at each sky position. Overlaid in the figure are the 
coverages of the SDSS DR7 data and the WISE preliminary release data. The overlap area 
between DR7 and the WISE preliminary release is 2344 deg 2 and defines the area analyzed 
in the rest of this paper. This area is slightly less than 30% of the total SDSS areal coverage. 

One characteristic feature of the WISE mission is that it does not have uniform sky 
coverage. For both the preliminary and the all-sky data release, the median depth-of-coverage 
is 15.65, 15.55, 14.85 and 14.84 at Wl, W2, W3 and W4 respectively (single exposure has 
11 seconds). 95% of the sky has coverage > 10.82 at Wl, and > 9.90 at Taking 10 

exposures per position as the depth-of-coverage in Wl and W2, and 9 in W3 and W4, the 
5a sensitivities are 17.3 mag (0.037mJy), 15.84 mag (0.079 mJy), 11.59 mag (0.67mJy) and 
8.0 mag (5.1 mJy) for the four band, respectively. 



1 See Section 6 in the Explanatory Supplement to the WISE All-Sky Data Release Products, 
http://wise2.ipac.caltech.edu/docs/release/allsky/expsup/index.html 
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Fig. 1. — Sky coverage of WISE and SDSS DR7 in Galactic Coordinates. The background 
map shows the WISE all-sky depth with colors indicating the number of single-exposure 
frames at each sky position. The white overlaid area shows the sky region covered by the 
WISE preliminary public data release, covering 57% of the full sky. The dark green region, 
primarily towards the North Galactic cap, is covered by the SDSS DR7 data. This paper 
analyzes data drawn from the overlap regions between the WISE preliminary public data 
release and SDSS DR7, a region covering 2344 deg 2 . 



2.2. WISE data 



We start with the WISE preliminary public release source catalog (Cutri et al. 2011) 
This catalog includes all sources detected in at least one band with photometric signal-to- 
noise ratio (SNR) > 7 after removing all artifacts. As described in the explanatory sup- 
plement to the preliminary data release products, source detection and photometric mea- 
surements are performed on all four bands simultaneously. The released catalog contains 
photometric measurements (or limits) in all four bands. 



2 A detailed description of the catalog can be found in the Explanatory Supplement to the WISE Prelim- 
inary Data Release Products at http://wise2.ipac.caltech.edu/docs/release/prelim/expsup/ . 
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WISE is the most sensitive at 3.4 /im and the least sensitive at 22 /im; the minimum 5cx 
sensitivities for these two bands are 0.05 and 6 mJy, respectively. For /„ oc z/*, the sensitivity 
ratio /„(3.4//m)// i ,(22/im) = [u (3 A jjim) / u (22 jim)} a = 0.0083 corresponds to a mid-infrared 
spectral index a ~ —2.56. Among all WISE VTl-detected galaxies at the limiting depth of 
that band, only 2% sources are also detected in W4, implying extremely red mid-infrared 
spectral energy distributions (SEDs), redder than /„ oc A 2,56 . 

The public released catalog does contain a very small fraction of sources which are 



detected only in W3 and/or Wi, but not in Wl and/or W2 (e.g., |Eisenhardt et al.||2012[ ). 
Due to the large sensitivity differences between Wl and W4, the entire catalog is essentially 
a 3.4 /zm-selected sample with SNR> 7 (<r wi ~ 0.15); only 0.16% of the catalog is very faint 
in Wl but bright in one of the other three WISE bands (e.g., SNR W i < 7, but SNR> 7 
for any of the other three bands). Therefore, we call this original catalog from the WISE 
data archive, without any additional cuts, the Wl sample. Many sources in this preliminary 
public release catalog have SNRjyi > 7, but low SNR in the other three WISE bands. When 
we study the color distributions, it becomes necessary to require the color errors to be small 
relative to the range of colors. Specifically, the Wl — W2 color for WISE galaxies roughly 
spans a narrow range, 0-2. We require the errors in Wl — W2 to be less than 0.27 mag, 
meaning we need a source catalog with SNRn/2 > 5 (corresponding to ay/2 < 0.22 mag). A 
similar argument is applied to the W2 — W3 colors, which have a larger spread of 0.5 — 5. 
In this case, we require SNRvi/3 > 3 (cr W3 < 0.35 mag). In summary, to obtain meaningful 
color distributions, we base our analyses below on the following three source samples: the 
original Wl sample (SNR^i > 7), a Wl/2 sample (SNRjyi > 7 and SNRw2 > 5), and a 
Wl/2/3 sample (SNR m > 7, SNRvk 2 > 5, SNR W3 > 3). 

Figure [2] shows the source fraction in each of the above three samples relative to the 
four-band merged, public released catalog. When applying additional photometric SNR cuts, 
the source sample sizes become increasingly smaller. At SNR > 5, 4.6 /im sources are only 
57% of the Wl sample. At SNR > 3, this fraction is only 15% and 2% at 12 /im and 22 [im, 
respectively. In the sense of source surface density, the Wl sample has ~ 8230 sources per 
deg 2 , whereas the Wl/2 and Wl/2/3 samples have only 4700 and 1235 sources per deg 2 , 
respectively. For comparison, the source surface density in the SDSS photometric catalog is 
about 27,300 sources per deg 2 , over three times higher than that of WISE 3.4 /zm. This is 
because SDSS detects many more low- luminosity (low-mass) galaxies than WISE (see, e.g., 



Donoso et al. 2012). 



For the three catalogs specifically selected for this paper, SNR^yi = 7 corresponds to ap- 
parent magnitudes between 16.8 and 17.6 mag (0.058- 0.028 mJy), SNRvk2 = 5 to 15.9 mag to 
16.5 mag (0.075- 0.043 mJy), and SNRn/3 = 3 to the limiting magnitudes of 12.2 to 12.8 mag 
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(0.38-0.22 mJy). 



WISE cannot measure absolute fluxes because the telescope down-links only the de- 
tector ramp slopes rather than each ramp readout. Therefore, WISE produces photometry 
relative to calibration stars, using the Vega magnitude system. The conversion factors to 
the AB system (m AB = m Ycga + A™) are 2.683 (Wl), 3.319 (W2), 5.242 (W3) and 6.604 
(W4); equivalently, zero magnitude corresponds to 309.5, 171.79, 31.676 and 8.36 Jy for 



the four bands, respectively (Wright et al. 2010) and (Jarrett et al. 2011). SDSS photom- 
etry directly output from the SDSS archive is in the asinh magnitude system. We trans- 
late asinh magnitudes into the AB system using information provided by SDSS^J When 
calculating SDSS— WISE colors, we convert SDSS photometry to the Vega system using 
r(AB) — r(Vega) = 0.16 (Fukugita, Shimasaku & Ichikawa 1995). Throughout this paper, 



WISE magnitudes are in the Vega system, SDSS magnitudes are in the AB system, but all 
colors are in the Vega system. 



2.3. SDSS data 



This paper uses the SDSS DR7 data (Abazajian et al. 2009), including the photomet- 



ric catalog, the main galaxy spectroscopic catalog (Brinchmann et al. 2004), the lumi- 



nous red galaxy (LRG) spectroscopic sample (Eisenstein et al. 2001) and the QSO sample 
(Schneider et al. 2010). The DR7 legacy survey catalog covers 8423 deg 2 and contains 



230 million sources. The main galaxy spectroscopic catalog is generated jointly by the 
Max-Planck-Institut fur Astrophysik and the Johns Hopkins University (MPA-JHU DR7 
main galaxy catalog; Brinchmann et al. 2004). It consists of almost 10 6 galaxies with 
Petrosian (1976) magnitudes brighter than r = 17.77 for which various derived physical pa- 



rameters are readily available. We adopt the following classification criteria as in Kauffman 
et al. (2003): star-forming (SF) galaxies are defined as having log[0 III/H/3] <1.3 + 0.61 



log[N II/Hq - 0.05] (Eq. 1 in Kauffmann etaL] 2003 ) ; AGNs (Seyfert and LINERs) are de- 



fined as having log([0 III]/H/3) > 1.19 + 0.61/[log([N II]/Ha) - 0.47] (Eq. 5 in Kewley et al. 



2001); and composite systems are defined as having log([0 III] /H/3) between the two values 
described by the above equations. These SF galaxy and AGN definitions are used in § 3.2 



when we utilize SDSS spectra to classify WISE sources. When based on WISE data alone, 
we do not use these definitions. Two LRG samples are selected in color-magnitude space, 
one to r = 19.2 (roughly volume-limited, to z — 0.38) and one to r = 19.5 (flux-limited, to 



3 See 
4 See 



http://www.sdss.org/DRl/algorithms/Iluxcal.html 



http:/ / www. mpa-garching. mpg. de/SDSS/DR 7/ 
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Relative Source Fraction to the WISE Public Catalog 




Band 

Fig. 2. — Source densities of our SNR selected WISE samples for the four bands at > 20°. 
The left-most bar, essentially equivalent to the 3.4 /zm SNR >7 catalog, has a source density 
of 8230 per deg 2 . The bars illustrate the source fractions with the corresponding SNR cuts 
at the four bands, relative to the original WISE preliminary public release source catalog, 
which includes all sources with SNR > 7 in at least one band. 



z = 0.55) (Eisenstein et al. 2001). The QSO sample is selected by their non-stellar colors 
or FIRST (Becker et al. 1995) radio emission to % — 19.1 for z < 3 and to i — 20.2 for 



3 < z < 5.5 (Richards et al. 2002). 



In addition, all SDSS photometric redshifts (when spectroscopic redshifts are not avail- 



able) used in sections below are the ones derived based on the neural network method (Oyaizu 



et al. 2008). We caution that, in general, broad band photometric redshifts are highly un- 



certain for strong type-1 AGNs/QSOs ( Assef et al.||2010 Brodwin et al.||2006 ). These strong 
type-1 AGNs from SDSS have the most complete spectroscopic redshifts. 



3. Analysis and Results 



3.1. Apparent magnitude and redshift distributions 

The WISE 3.4 and 4.6 /im bands primarily sample emission from stellar photospheres, 
whereas the 12 and 22 /im bands are more sensitive to dust emission heated by stars and 
accreting black holes. Therefore, WISE readily identifies a variety of galaxy populations 
in the mid-infrared sky. To characterize these populations, we first address several basic 
properties of the matched W/SE-SDSS catalog, including optical brightness distributions, 
redshifts, the fraction of WISE sources undetected or with very faint optical counterparts, 
and optical thru mid-infrared colors of matched sources. 

We carried out source matching between WISE and SDSS DR7 with a matching radius 
of 3". This is based on the fact that the WISE angular resolutions are 6'.'1, 6'.'4, 6'.' 5, and 12" 
for the four bands, respectively. Figure [3] shows the percentage of WISE sources with SDSS 
r-band optical counterparts for the four WISE bands. To a depth of r = 22.6, corresponding 



to the 50% completeness depth of SDSS (Abazajian et al. 2004), 91% of the Wl/2/3 and 



96% of the Wl/2/3/4 (SNR^ 4 > 3) samples have optical counterparts. In contrast, only 
72% and 86% of the Wl and the Wl/2 samples have optical counterparts with r < 22.6. 
This implies that the bulk of 12 and 22 /im galaxies are dusty star-forming galaxies and 
AGNs at low redshifts [z ~ 0.1 — 0.3), whereas 15-25% of 3.4 and 4.6 /im sources could be 
massive early-type galaxies which are fainter than the SDSS photometric limits and are at 



z>l. Indeed, Gettings et al. (2012) reports on a WISE-selected galaxy cluster at z = 0.99 



which is well-detected in Wl and W2 but is undetected by SDSS. 

Figure [4] shows the optical brightness and photometric redshift distributions for the 
mid-infrared sources detected in the four bands. It is clear that WISE 22 /im sources have 
corresponding optical magnitudes which peak at r ~19, one magnitude brighter than the 
peaks for the Wl, W2 and W3 samples. For comparison, the parent SDSS photometric 
catalog continues to rise until completeness causes an apparent drop at r ~ 22. This confirms 
that SDSS detects many more low-mass, low-luminosity galaxies than WISE and the bulk 
of WISE galaxies with optical counterparts are relatively local, at 0.1 < z < 0.3. 

About 28% and 14% of WISE 3.4 /im and 4.6 /im sources are very faint or without 
any optical counterparts in the SDSS photometric catalog. These optically faint 3.4 /im 
galaxies are predominantly massive early-type galaxies at z > 1, though heavily obscured 
galaxies and AGN also contribute at some level. Figure [5] examines the Wl — W2 colors of 
this optically faint population, showing that Wl sources with very faint r-band magnitudes 
(r > 22.6) have redder Wl — W2 colors than optically brighter Wl sources. As discussed 



below and shown with color-redshift tracks for galaxy templates from Asscf et al. (2010), 
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redder Wl — W2 color is an indicator of early-type and star-forming galaxies being at higher 
redshifts, z > 1. 
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Fig. 3. — Percentage of WISE sources matched with SDSS photometric catalog as a function 
of WISE magnitude. 

What types of galaxies are detected by WISE at 3.4 and 22 microns? One way to 
answer this is to calculate expected WISE magnitudes as a function of redshift using em- 
pirical galaxy templates. Figure [6] shows the expected Wl and W2 magnitudes for an 
early-type galaxy template, as well as W3 and W4 magnitudes for an infrared-bright galaxy 
(IRAS 19254—7245) as a function of redshift. This figure suggests that ignoring evolution, 
L* early-type galaxies should be visible in WISE 3.4 /xm data out to redshifts of 1.5-2. An- 
other important result emphasized by Figure [6] is that the observed Wl magnitudes do not 
change significantly over z ~ 0.5 — 1.5 for early-type galaxies; this benevolent /c-correction 
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Fig. 4. — The r-band magnitude and photometric redshift distributions for the WISE-SDSS 
matched sources from the four WISE samples. For both panels, the distributions for Wl 
and Wl/2 samples are indistinguishable, with the Wl sample in blue and the Wl/2 sample 
is green. In the left panel, for the Wl/2/3/4 sample, we scaled the y-axis by a factor 10 
for visual clarity. Similarly, in the right panel, the y-axis for the Wl/2/3/4 sample is scaled 
by 2. For comparison, we also plot the r-band magnitude distribution for the full SDSS 
photometric catalog (scaled down by a factor of 1000). Note that SDSS optical magnitudes 
are in the AB system, calculated from SDSS asinh magnitudes. The r-band magnitude 
histograms have a bin size of 0.2 mag. The majority of 12 /mi sources have bright SDSS 
optical counterparts and are at low redshifts (z ~ 0.2). 



has been noted repeatedly for the similar 3.6 /mi band of Spitzer (e.g., Eisenhardt et al.||2008 



Mancone et al. 2010 Galametz et al. 2012). In contrast, optical r-band magnitudes steeply 



decline in brightness with redshift. 

Figure [7] illustrates how Wl and W4 magnitudes change with optical brightness, r, 
and photometric redshift. Wl traces r until r ~ 20 — 21, at which point Wl magnitudes 
remain within a narrow range of 15.8 - 16.8 mag, while optical brightness becomes increasingly 
fainter, stretching over three magnitudes. Figure [6] and Figure [7] illustrate the same effect, 
i.e., most galaxy SEDs have very steep slopes at ultraviolet through optical wavelengths, 
and much shallower slopes in the near-infrared. The turn-over at r ~ 20 in Figure [7] implies 
that the bulk of Wl sources are low-redshift (z < 0.5), typical SDSS-detected galaxies, and 
Wl sources with r > 20 tend to be at higher redshifts (z ~ 0.5 — 2). 
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Fig. 5. — Normalized Wl — W2 color distribution for Wl/2 sample with bright optical 
counterparts (solid line; r < 22.6) compared to sources with faint or undetected optical 
counterparts in the SDSS data (r > 22.6). 



3.2. Separating powerful AGN/QSOs at z < 3 from galaxies 

With four-band mid-infrared photometry, can WISE colors alone provide diagnostics 
for different types of sources as well as crude redshift information? Figure [8] addresses this 
question by showing Wl — W2 vs. W2 — W3 colors for the Wl/2/3 sample. The three high 
source density regions are: (1) the stellar locus with both colors near zero; at mid-infrared 
wavelengths, emission from Galactic stars from spectral class O through ~ T3 is dominated 
by the Rayleigh- Jeans tail of the blackbody spectrum, thereby yielding Vega-system colors 
near zero for most Galactic stars (note that pure elliptical galaxies with no dust at z < 0.1 
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Fig. 6. — Predicted r-, W1-, W2-, W3- and iy4-band magnitude as a function of redshift 
(corresponding to black, blue, green, magenta and red lines respectively) using two empirical 
galaxy SED templates, elliptical galaxy and IRAS19254-7245. For the early type galaxy, the 
template is normalized to an L* galaxy at z = 0.043 with Lk = 3.5 x 1O 11 L (Lin et al.||2004 



Mancone et al. 2010). Minimum WISE sensitivity limits are shown as dashed horizontal 



lines. This shows that WISE 3.4 /xm images are sensitive to typical early-type galaxies out 
to z ~ 1; this sensitivity is enhanced for more luminous early- type galaxies and/or higher 
ecliptic latitude fields where the WISE coverage is deeper. 



also sit within the stellar locus); (2) the red Wl - W2 cloud with Wl - W2 > 0.8; and 
(3) the bluer Wl — W2 sequence, spanning a wide range in W2 — W3 color. Here we do 
not see any extremely cool brown dwarfs with very red Wl — W2 color due to methane 
absorption since such sources are very faint in W3, and thus are excluded in our Wl/2/3 



sample (Kirkpatrick et al. 2011). 



What is the physical basis of Figure ^ Figure [9] through 11 examine this question from 
different angles. Figure [9] shows empirical galaxy and AGN templates from 0.03 — 30 fim 
based on ~ 20000 objects with multi-wavelength data in NOAO Deep Wide-Field Survey 
(Assef et al. 2010). For systems with strong nuclear heating, the mid-infrared SEDs tend 
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Fig. 7. — WISE 3.4 /im and 22 /zm magnitude as a function of optical r-band magnitude. The 
colors indicate photometric redshifts taken from the SDSS catalog. As expected, optically 
fainter galaxies, i.e., sources with redder r — Wl or r — W4 colors tend to be at higher 
redshifts. 



to be roughly a rising power-law, resulting in Wl — W2 > 0.8. For systems with dominant 
stellar emission at z < 1, this color is smaller because the 3.4 jum band is sampling below the 
1.6 //m peak of the stellar emission, making the Wl — W2 color bluer; at z > 1, the 3.4 /zm 
band samples below the near-infrared stellar peak, getting fainter, thereby providing redder 
Wl - W2 colors. 



Figure 10 shows the Wl — W2 and W2 — W3 color tracks as a function of redshift 
calculated using the set of SED templates of Assef et al. (2010), the Arp220 SED template 



of Pollettaet al. (2007) and the IRAS15250+3609 template of Vega et al. (2008). The model 



colors clearly do not span the full range of observed values in Figure [8] The assumed SED 
templates, generated from optical and Spitzer observations, are by no means complete. These 
templates include elliptical galaxies, star forming Sbc and Im type galaxies, type-1 AGNs 
and local ULIRGs Arp220 and IRAS15250+3609. The AGN tracks have three different dust 
obscuration factors; dust-obscured ones are type-2 AGNs. 



As opposed to Figure 10 which shows WISE colors based on SED templates, Figure 11 is 
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the color-color plot for known sources. Specifically, we plot SDSS spectroscopically classified 
sources, including star-forming (SF) galaxies, galaxies hosting AGN, SF/AGN composite 
systems, LRGs, QSOs and z <: 2 WISE-selected ULIRGs. These figures enable us to conclude 
that WISE colors alone are effective in separating strong AGNs from star-forming galaxies, 
and in separating stars from extragalactic sources (other than early-type galaxies at very 
low-redshift). 



Figure [8] through 11 all suggest one unique application of WISE data over the entire 
sky — selecting strong AGN/QSO candidates at z < 3 using Wl — W2 color. Using mid-IR 
colors to select AGNs has been noted earlier in Spitzer studies over relatively small areas and 



with much smaller AGN samples (e 
W1 — W2 method, briefly discussed 


.g., Lacy et al. 


2004 


Stern et al. 


2005). The selection of 


in Jarrett et al. 


(2011 


), is tested and discussed in detail in 


Stern et al. 


(2012) and Assef et al. ('. 


2012 ) (see also i 


Vshby et al.|2009 


Assef et al.|2010; Eckart 


et al.|2010 


Edelson & Malkan||2012 


Massaro et al.||2012 


Wu et al.||2012 Mateos et al.|2012). 



With the current analysis, we, for the first time, investigate this simple selection criterion 
using the wide-area, large spectroscopic database of SDSS. Figure 12 plots Wl — W2 vs. 
W2 — W3 for SDSS spectroscopically confirmed QSOs, color-coded by redshift. Normal, star- 
forming galaxies are also plotted. We see that QSOs at z < 3 mostly have Wl — W2 > 0.8, 
and are reasonably well separated from star-forming galaxies. However, z ~ 3 — 5 QSOs 
have bluer Wl — W2 colors due to the 3.4 and 4.6 /xm filters sampling rest-frame optical 
wavelengths, making such quasars difficult to distinguish from normal galaxies based on 
mid-IR colors alone. A similar issue with Spitzer selection of distant quasars was pointed 



out by Assef et al. (2010). The selected AGN candidate sample has two unique features: (1) 
it provides a more complete selection of QSOs, particularly at z ~ 2 — 3, where the SDSS 



optical color method has issues (Richards et al. 2002 Schneider et al. 2010); and (2) it is 
less affected by dust than optical selection methods, and thus is more sensitive to obscured, 
type-2 AGNs. 

In addition to the color criterion of Wl — W2 > 0.8, another important condition 
for selecting AGN with high completeness and reliability is W2 magnitude and SNR limits 
(see Stern et al. (2012) and Assef et al. (2012) for the full details). Using the deep, multi- 



wavelength data in Cosmic Evolution Survey (COSMOS; Scoville et al. 2007), Stern et 
al. (2012) found that Wl - W2 > 0.8 in combination with W2 < 15.05 at SNR > 10 
identifies 78% of Spitzer mid-IR AGN candidates (as the truth sample) with a reliability of 



o. 



Assef et al. (2012) uses higher latitude WISE data in the wider-area Bootes field to 
further investigate the completeness and reliability of the WISE AGN selection as functions 
of magnitude and SNR. As shown in Figure [10] and discussed in § 3.1, optically faint, high- 



redshift elliptical and Sbc galaxies tend to have very red Wl — W2 colors. These sources 
are contaminants to WISE AGNs based on a pure Wl — W2 > 0.8 color cut. To limit the 
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contamination, Stern et al. (in prep) and Assef et al. (2012) found that W2 < 15.2 at 
the number of sky coverage > 15 can achieve completeness of 74% and reliability of 89%. 
Here similar to Stern et al. (2012), the reliability and completeness of WISE AGN sample is 
derived based on using the Spitzer color selected mid-IR AGNs as the truth sample. For our 
chosen WISE + SDSS overlap sub-region, more than 80% of the sky have the coverage > 15, 
and 99% have the coverage > 10. At W2 = 15.2, the median source SNR and the number sky 
coverage are 12 and 22 respectively. Therefore, for our analysis of WISE AGNs throughout 
this paper, we adopt W2 < 15.2 criterion. Our AGN selection criteria of Wl — W2 > 0.8, 
W2 < 15.2 should maintain a reasonable completeness and a good reliability. 

What is the source density of WISE-selected sources? In the large, high Galactic latitude 
WTST^-SDSS overlap area we analyze here, the Wl/2/3 sample contains 12% Galactic stars, 
12% luminous AGN candidates, and 70% normal galaxies (possibly containing weak AGNs). 
The remaining 6% of W3 sources are extremely red, with either W2 — W3 > 4.5 or Wl — 
W2 > 1.8. These present a rare population of red objects, and are likely dusty systems at 
high redshift (see 



3.4). With the selection of Wl 



W2 > 0.8, W2 < 15.2, the WISE 
AGN/QSO candidates have a surface density of 67.5 ± 0.14 per deg 2 . Considering only the 
optically bright subsample, with r < 21, the AGN candidate surface density is 35 per deg 2 . 
For comparison, Stern et al. (2012) finds 61.9±5.4 per deg~ 2 AGN candidates with 95% 
reliability for W2 < 15.05, and Assef et al. (in prep) finds 137±4 per deg -2 candidates for 
W2 < 17.11. 

SDSS spectroscopic QSOs are optically bright sources. What is the fraction of this 



population detected by WISE 3.4 and 22 /zm? Figure 13 addresses this question. Overall 



WISE detects a high fraction of SDSS QSOs from DR7 Schneider et al. (2010) catalog, and 



z > 5.7 QSOs from Fan et al. (2006); Jiang et al. (2009), Willott et al. (2009) and Willott 



et al.| ( |2010[ ). Of the entire QSO sample, only 18.9% have WA detections at SNR> 5a, and 
much higher fraction, 89.8% have Wl detections at SNR> 7a. We note that WISE Wl 
is quite sensitive in detecting z > 5.8 QSOs, with a detection rate of ~50%, as shown in 



the bottom panel of Figure 13 In particular, the highest redshift quasar currently known, 
ULAS J1120+0461 at z = 7.085 (|Mortlock et^[20lT|, is also detected by WISE 3.4 /mi. 



Blain et al. (2012) presents a detailed discussion of the WISE properties of z > 6 optically 



selected QSOs. 



3.3. Type-2 AGN candidates 



WISE is very sensitive to dust-obscured objects. This naturally leads to two questions: 
(1) can WISE colors be used to identify type-2 AGNs; and (2) what is the relative fraction of 
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Fig. 8. — Color-color distribution of WISE sources from the 12 /xm-selected Wl/2/3 sample. 
To illustrate the large range in source density in color-color space, we combine a grey-scale 
plot for high concentration regions (see scale-bar on right) and individual points for low 
source density regions. 



type-1 and type-2 QSO/AGNs selected by WISE? Based on actual data, Figure 14 illustrates 
that type-1 and type-2 QSO/AGNs likely have different mid-infrared color distributions by 
comparing the Wl — W2 vs. W2 — W3 colors of unresolved SDSS sources (point sources, 
including QSOs) with extended sources. To ensure reliable star/galaxy classification, we 
require r < 21 (e.g., 95% reliability quoted in Table 2 in |Stoughton et al.|2002 ). WISE colors 
of unresolved sources (e.g., Galactic stars and quasars) are clearly separated from extended 
galaxies. Furthermore, in terms of unresolved SDSS sources, note that the Galactic stars (at 
the origin) are clearly distinct from the unresolved quasars (with redder mid-infrared colors) . 



One cautionary note about type-2 AGN is that any discussion and conclusions will 
depend how such sources are defined. In the current analysis, we define two type-2 AGN 
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Fig. 9. — Empirical templates for quasars and normal galaxies from Assef et al. (2010). The 
major difference between galaxies and AGN is at 1-2 /mi: AGN have a minimum in that 
wavelength range, while galaxies have a peak in that wavelength range. Star-forming and 
passive galaxies are easily distinguished from their longer- wavelength data (> 4 /im). 



samples. The first are AGN candidates with WI-W2 > 0.8, W2 < 15.2, and extended SDSS 
r-band morphologies (SDSS TYPE = 3; r < 21). Unfortunately, this definition only works 
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for optically bright AGN candidates with spatially resolved host galaxies in the SDSS data. 
It is also worth noting that mid-IR color criteria for AGNs tend to select systems with strong 
black hole accretions. For fainter type-2 AGN candidates, we use an alternative definition: 
Wl — W2 > 0.8, W2 < 15.2 and very red optical-to-mid-infrared color — i.e., optically 
faint, WISE 4.6 /im AGN candidates. As emphasized earlier, when selecting WISE AGNs 
using Wl — W2 colors, the W2 < 15.2 criterion is very important for limiting contamination 



from high-redshift galaxies to less than 20% (Assef et al. 2012). 



For the former, optically bright type-2 AGN definition, we estimate ~ 31% of r < 21 
WISE-selected AGN/QSOs to be spatially resolved. This suggests an almost 2:1 type-1 to 
type-2 AGN ratio. Since type-2 AGN here only include bright AGNs with well-resolved 
host galaxies in optical images, it is not surprising that the type-2 to type-1 ratio is so 
e.g., lower than what is required to explain the hard X-ray background (Gilli et al. 



low 



20071 |Comastri et al.||2011[ ). With SDSS + WISE data, it is also possible to use r - W2 
color to select type-2 AGNs. Figure 15 examines this technique, showing r — W2 colors 
for both type-1 and optically bright type-2 AGNs with extended SDSS morphologies. The 
morphological criterion imposes a strong selection function, producing an artificial deficit 



of luminous type-2 AGN in Figure [15j Higher luminosity sources will tend to be at higher 
redshifts and high-redshift obscured AGNs drop below the r = 21 limit imposed to provide 
reliable morphologies. Here L 12Atm is calculated using spectroscopic or photometric redshifts, 
including small fc-corrections (less than 0.1 magnitude) based on the appropriate type-1 or 



type-2 AGN template from Assef et al. (2010). Type-1 and type-2 AGNs clearly have very 
different r — W2 color distributions, especially at the high luminosity end. This is consistent 



with what has been found by Hickox et al. (2011) based on X-ray and Spitzer data in the 
Bootes field. This result lends support for using r — W2 to select type-2 AGNs which do not 
have reliable optical morphology information. 



Furthermore, the feasibility of this method is illustrated in Figure [16j showing the 
calculated optical/mid-IR colors based on assumed sets of SED templates as a function of 



redshift. These galaxy templates are constrained by actual Spitzer observations (Polletta et 



al.|2007 Vega et al.|2008 ). We see a strong divergence in optical/mid-IR colors of type-1 and 
type-2 QSO/AGNs at z > 0.5. Figure 17 further examines this technique, showing r — W2 
vs. Wl — W2 for WISE sources with W2 < 15.2, matched with SDSS r-band data down to 
the faintest limits. The secondary branch in the figure shows that sources with red r — W2 
colors are potential type-2 AGNs. If we use the criteria of Wl - W2 > 0.8, W2 < 15.2, 
and r — W2 > 6 to define a sample of red AGN as potential type-2 sources, we find that 
type-2 candidates account for ~ 29% of all WTS'.E-selected AGN candidates to that depth. 
The surface density of such a type-2 AGN candidates is about 16.4 per deg 2 , whereas r < 21 
bright with resolved morphology, type-2 AGN is about 31% of all AGNs, yielding surface 
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density of 17.3 per deg . Stern et al. (2012), using the Hubble Space Telescope imaging 



available in the COSMOS field, find that ~ 50% of WISE-selected AGN candidates to a 
depth of W2 = 15.2 are unresolved in / 8 i 4 . 

Our intention here is to illustrate the methods of using colors to select large samples 
of type-2 AGN candidates over wide areas of sky. Our crude estimate of the ype-2 AGN 
fraction (integrated) is on the order of l/3rd of all AGNs. However, we caution that to 
truly understand the implications of these numbers, we will need to have the redshift and 
luminosity information. Such studies utilizing WISE plus other ancillary spectroscopic data 



in smaller regions of sky are included in Assef et al. (2012) and Stern et al. (in prep) 



3.4. High-redshift ULIRG candidates 



One of the primary WISE mission goals is to identify extremely luminous (Ljr > 
10 12 ~ 13 L & ) dusty starbursts and AGNs (i.e., ULIRGs and hyper-luminous infrared galaxies, 
or HyLIRGs) at high redshift. We have pursued two basic approaches in finding high-redshift 
ULIRGs. One is to use WISE data alone, and the second is to combine WISE with optical 
data. 

Using WISE colors alone, high-redshift ULIRGs can be identified by requiring very 
red, rising mid-infrared SEDs, redder than the spectral index a ~ —2.56 corresponding 
to the WISE sensitivity limits (see § 2.2). Specifically, this means significant detections 



(SNR > 3 - 5) in W3 or W4, but no detections in Wl and W2, i.e., so-called PUlIU2-dropouts 



(Eisenhardt et al. 2012). Follow-up optical spectroscopy of more than 100 candidates using 



the Keck, Gemini, Magellan, and Palomar telescopes have demonstrated that the majority 
of these candidates are indeed at redshifts of 1-3, implying very high luminosities given the 



high fluxes at 12 and 22 /im (Eisenhardt et al. 2012 Wu et al. 2012). This type of ULIRG 



is rare, with a surface density of just 0.02 sources per deg 2 . One interesting feature of these 
ULIRGs is their extremely steep mid-IR spectral slopes. Comparing Figure 10 and 11, we 
see that the newly discovered WISE-selected z ~ 2 ULIRGs have much redder W2 — W3 
colors than the calculated high- z colors using local ULIRG SED templates, including the 



reddest one such as IRAS15250+3609. Eisenhardt et al. (2012), Wu et al. (2012), Bridge et 



al. (2012) and Yan et al. (2012) present complete discussions on the PUlIU2-drop selection, 



follow-up spectroscopy, as well as far-infrared photometry from CSO and Herschel for some 
of these extremely luminous {Lm ~ 1O 13_14 L ) ULIRGs at z ~ 2. 



The second method is to utilize optical/mid-infrared colors. As shown in Figure 16 
r — Wl color becomes redder at higher redshift. This is true as well for the other three WISE 
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bands. Indeed, there is a rich literature using red optical-to-mid-infrared colors with Spitzer 
to select high-redshift galaxies. For example, red 24 /im to r-band colors (r — [24] > 14; 
Vega) have been used in the Extragalactic First Look Survey (XFLS; (Yan et al. 2004)) 
and Bootes field to successfully select many highly obscured galaxies at z ~ 2 (Yan et al. 



2005 


2007 


Dey et al. 


2008) 



such as IRAS 19254-7245, Mrk 231, UGC5101 and IRAS 08572+3915 have very red colors 
at any redshifts. Using r — W4 > 14 as the selection criterion, WISE identifies 0.9 ± 0.07 
high-redshift ULIRG candidates per deg 2 for SNRyi/4 > 5. This signal-to-noise ratio cut 
roughly corresponds to a 22 fim flux density > 2.5 mJy This is shown in Figure 18 If we 
instead apply a uniform flux density cut of f u (22fim) > 5 mJy, the surface density of ULIRG 
candidates is 0.41 ± 0.05 per deg 2 , comparable to that of Dey et al. (2008) using Spitzer 
data over the ~10 deg~ 2 Bootes field with the flux density cut of /„(24//m) > 5 mJy. To 
verify our analysis, we chose a subset of the SDSS-WISE overlap region covering 180 sq. 
degrees. All of the red, high-redshift candidates with r — WA > 14 were visually examined. 
We found a small percentage of contaminants (8%). Our derived surface density values have 
been corrected for this percentage. The WISE all-sky data offers an excellent opportunity to 
assemble a large sample of 22 /zm-selected "dust-obscured galaxies" , or DOGs, at z > 1 — 2. 



3.5. Keck spectroscopy of a sample of high-redshift ULIRG candidates with 

very red r — WA color 

The WISE team has carried out optical spectroscopic follow-up observations of high- 
redshift ULIRG candidates. These candidates include both WlIV2-dropouts as well as 



extremely red optical/mid-IR sources. Eisenhardt et al. (2012) focused on one particular 



IVlW2-dropout, J181417.29+341224.9 at z = 2.452, whose optical spectrum shows a typical 
star-forming galaxy whereas its IR SED suggests a highly obscured AGN with a bolometric 
luminosity of 3.7 x 1O 13 L . Wu et al. (2012) presents optical spectra of a sub-set of W1W2- 



dropouts with millimeter observations, making a comparison between IVlIV2-dropouts and 
Spitzer-discovered dust obscured galaxies ("DOGs"). 

The follow-up observations included a sample of WISE sources with r — W4 > 14. 
Table [T] lists the source information for eight galaxies with good spectra. These spectra were 
taken using Low Resolution Imager and Spectrograph (LRIS,( Oke et al.||1995 )) on the Keck 
telescope during the nights of March 10, April 10, and May 10, 2011. The total on-target 
integration times range from 10 to 15 minutes. Figure [19] shows the wavelength calibrated 
spectra in the observed frame for these 8 sources. Prominent emission lines include Lya, 
C IV A1549, [ionMg2] A2800, Ne V A3426, [O II] A3727, and [O III] A5007 etc, strong 
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spectral features. W1422+5613 has a C IV A1549 FWHM of 1480 km/s, not accounting 
for the spectral resolution. If we adopt the cutoff at 2000 km/s for the type-1 and type-2 
separation based on C IV A1549, this object is indeed a type-2 AGN. This cutoff value is 
very conservative comparing to the mean C IV A1549 velocity width of 5600 km/s, estimated 



the SDSS QSO sample flShen et al.||2008[ ). Similarly, W145 7+2932 is a type-1 AGN with a 
C IV A1549 FWHM of 3940 km/s. W1603+3627 has a C IV A1549 FWHM of 1790 km/s, 
however, its C III] A 1909 FWHM is broad, 3170 km/s and the [ionMg2] A2800 is broad as 
well. Since C IV A1549 line is a complex line and could have some absorption. We classify 
this object as type-1 based on C III] A 1909 and [ionMg2] A2800 instead. W1604+2041 is a 
type-2 AGN with a Ne V A3426 FWHM of 1350 km/s. W1617+3355 has a ratio of [O III] 
A5007 to H/3 of 3.7, which is in a range that could be either AGN or star-forming based on 
the Baldwin, Phillips & Terlevich (BPT) diagram. Its [O III] A5007 FWHM is 610 km/s, 
making it consistent with either a type-2 AGN or a star-forming galaxies. W1617+2336 is 
a type-1 AGN with a C IV A1549 FWHM of 4800 km/s. W1636+3309 has a ratio of [O III] 
A5007 to H/3 >20, which implies an AGN, however, the available spectral line widths for 
a small number of features are narrow, suggesting it being a type-2. The classification for 
this source is not clear. W1708+3315 has [O III] /H/? = 3.0, again making it ambiguous, 
either a type-2 AGN or a star-forming galaxy. Table [T] lists these classifications. Although 
the number of spectra is small and there are some ambiguous systems, the source types 
among the randomly selected 8 targets are dominated by type-2 AGNs. This supports the 
interpretation that the high mid-IR fluxes at 22 /im are due to dust emission heated by 
obscured central AGNs. 

We plot the colors and brightness of these 8 sources as a function of their spectroscopic 



redshifts. The three panels in Figure 20 show tends which suggest that optically fainter or 
r — W4 redder candidates could be at higher redshifts. The brightness in the 22 /im band 
does not appear to correlate with redshift, with many high-redshift ULIRGs being quite 
bright in W4. 



4. Discussion and Summary 

With the public data releases, WISE has delivered to the community the most sensitive 
all-sky mid-infrared map of our generation. For extragalactic research, what can WISE 
data provide? In this paper, we present a phenomenological study that characterizes the 
observational properties of mid-infrared extragalactic sources and identifies color selection 
criteria for isolating large samples of QSOs, dust-obscured type-2 AGNs and luminous high- 
redshift ULIRG candidates. 
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With 5a sensitivities < 0.05, 0.1, 0.75 and 6mJy at 3.4, 4.6, 12 and 22 fan, the Wl, 
Wl/2, Wl/2/3 and Wl/2/3/4 samples have source surface densities of 8230, 4700, 1235 and 
150 per deg 2 , respectively. At the limit of the data, only very red mid-infrared sources with 
spectral slopes steeper than ~ —2.56 are simultaneously detected in both Wl and W4. The 
Wl source density is one-third that of the SDSS photometric catalog, suggesting that WISE 



3.4 /im images miss many low- mass, low- luminosity, blue galaxies (Donoso et al. 2012). We 
find that 28% of the Wl sample have r-band magnitudes fainter than 22.6, including many 
sources lacking optical counterparts. We present observational evidence that suggests that 
optically faint 3.4 /xm sources are likely early-type galaxies beyond the redshift limit of the 
SDSS imaging data. 

WISE 3.4 /zm data is sensitive to bright quasars. Of the entire SDSS optical QSO sample 



(Schneider et al. 2010), 89.8% have Wl detections at SNR> 7. Known to date, WISE even 



detects 3.4 /im emission from the highest redshift QSO, ULAS Jl 120+0461 at z = 7.085 



(Mortlock et al. 2011). In contrast, only 18.9% have WA detections at SNR> 5. Some of 
these optically selected bright QSOs detected at 22 /im are not just bright in all bands, but 
indeed have strong mid-IR dust emission. 

The unique advantage of WISE all-sky data for extragalactic sources is the diagnostic 
powerful of its mid-infrared colors. This power comes from the simple fact that pure stellar 
systems, star-forming galaxies and QSO/AGN have distinctly different near-infrared SEDs 
(see Figure [9J, thus yielding very different observed mid-infrared colors for different types of 
objects. We show that WISE colors alone can separate source populations, including Galac- 
tic stars, star-forming galaxies and QSO/AGN. We present three useful applications of WISE 
all-sky data: (1) We select QSO/strong AGN at z < 3 using W1 — W2 color. This population 
of QSOs is particularly interesting to inventory since a large fraction of QSOs at z ~ 2 — 3 are 
missed by the SDSS optical color selection. Beyond z ~ 3, the WISE QSO color selection 
starts to fail as the 3.4 /im and 4.6 /zm bands sample wavelengths blue-ward of lA, and the 
observed WISE colors start to have overlap with low-redshift star- forming galaxies. (2) We 
demonstrate the possibility for selecting type-2 AGN/QSO candidates using W1 — W2 > 0.8, 
W2 < 15.2 and r — W2 > 6 (or SDSS optically resolved morphologies for r < 21 sources). 
The WISE data allow a more complete census of dust-obscured, actively accreting, super- 
massive black hole systems than allowed by the current generation of deep X-ray surveys. 
(3) High-redshift [z > 2) dust-obscured, extremely luminous ULIRGs/HyLIRGs are another 
population of extragalactic objects which the WISE data can be used to identify. With 



follow-up optical spectroscopy, Eisenhardt et al. (2012) demonstrates the efficient selection 



of z ~ 2 — 3 ULIRGs with extremely red mid-infrared slopes (the W / lW / 2-drop-out method). 
In addition, we show that extremely red optical-to-mid-IR colors, e.g., using r — WA > 14, 
can also be used to select high-redshift dust-obscured starbursts and AGN. The candidate 
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surface density for this later selection is 0.9 ± 0.07 at SNR^ 4 > 5. This is consistent with 
the number inferred from the previous studies using Spitzer data over much smaller areas. 

Optical spectroscopic follow-up of a small number of r — WA >14 ULIRG candidates 
confirms that the color selection indeed works, identifying IR luminous galaxies over a wide 
range of redshift, z ~ 0.7 — 3. Despite small number statistics, we find indications that 
optically fainter and r — W4 redder sources tend to be at higher redshifts. The Keck optical 
spectra for the 8 sources detect many typical strong emission nebular lines seen among AGNs, 
such as Lyct, C IV A1549, Ne V A3426, [O III] A5007 lines. In addition, of the 8 sources 
with the optical spectra, five have spectral type of type-2 AGN, one type-1 AGN with fairly 
broad C IV A1549 emission line, and two possible star-forming galaxies. 
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Fig. 10. — Color-color tracks as a function of redshift for several galaxy and AGN templates 
taken from Polletta et al. (2007), Vega et al. (2008) and Assef et al. (2010). The blue 
tracks are for AGNs with target symbols indicating z — 0, black dots for z = 2,4 and open 
squares for z — 6. The AGN tracks are shown with three different dust obscuration factors. 
The red tracks are for the local ULIRGs Arp220 and IRAS15250+3609 with target symbols 
for z — 0, solid dots for z = 1,2,3 and open squares for z = 4. The magenta tracks are 
for normal galaxy templates with target symbols for z = 0, solid dots for z = 1 and open 
squares for z = 2. 
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Fig. 11. — As per Figure [8j with source classifications indicated, based on the SDSS DR7 
spectral measurements from the MPA-JHU catalog (see text for details). Spectroscopic 
confirmed QSOs and LRGs are also plotted. Additionally, we plot several extremely red 



WISE-selected ULIRGs at z > 2 from Eisenhardt et al. (2012), Stern et al. (in prep), Bridge 



et al. 


(2012 


)• 


Wu et al. 


(2012 


), and 


Yan et al. 


(2012 
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[4.6] — [12] (Vega mag) 

Fig. 12. — WISE color-color distribution for SDSS spectroscopically confirmed quasars, 
color-coded by redshift. Grey points are normal star-forming galaxies. A simple Wl — W2 
color cut efficiently and robustly separates quasars at z < 3 from normal galaxies. Higher 
redshift quasars, however, enter into the galaxy locus in this color-color diagram. 
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Fig. 13. — The top two panels show the WISE magnitudes (Wl and W4) as a function 
of redshift for the SDSS QSO sample (Schneider et al. 2010[ ). As expected, more QSOs 
are detected at 3.4 /im than at 22 /mi. The bottom panel shows the percentages of WISE 



detected SDSS QSOs as a function of z. The z > 5.7 sources are from Fan et al. (2006), Jiang 



et al. (2009) and CFHT surveys (Willott et al. 2009 2010). 
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Fig. 14.— WISE colors for SDSS unresolved sources (top panel; SDSS TYPE=3) and 
extended galaxies (bottom panel; SDSS TYPE=6). Comparison of the two panels shows 
that SDSS point sources have very different Wl — W2 and W2 — W3 color distributions 
from that of extended objects. In the top panel, the solid contours represent the distribution 
of the extended objects shown below. Furthermore, unresolved SDSS sources clearly separate 
into their two primary constituencies, stars at the origin and quasars with redder colors. 
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Fig. 15. — Rest-frame 12 /iiri luminosity for the optically bright (r < 21) AGN sample, 
including small k- correct ions calculated using appropriate SEDs for type-1 (blue) and type- 
2 (red) AGN. The differentiation between AGN type is based on optical morphology. The 
lack of high-luminosity type-2 AGN is due to the selection effect which uses SDSS extended 
morphologies as a criterion. More luminous type-2 AGN will be at higher redshifts where 
the host galaxies brightness falls below the optical magnitude limit imposed here to ensure 
reliable morphologies. Likewise, the host galaxies of low-luminosity, broad-lined (type-1) 
AGN are likely detected by SDSS, and thus such systems are classified as type-2 AGN using 
our morphological criterion. 
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Fig. 16. — Optical-to- WISE color as a function of redshift for a set of galaxy templates. 
Red ULIRGs refer to local, highly dust-obscured ULIRGs with red mid-infrared SEDs, e.g., 
Mrk231 (a type-1 AGN), UGC5101, IRAS 0872+3915 and IRAS 19254-7245, taken from 



Polletta et al. 


( 


2007) and Vega et al. 


( 


2008) 



-33- 




Fig. 17. — Color-color diagram, plotting r — W2 vs. Wl — W2. Vertical dashed lines 
show cuts at various Wl — W2 color; the plot on the right shows the corresponding relative 
r — W2 color distribution for sources redder than those cuts. The right-hand side panel is 
the histogram of the r — W2 color distribution; and the x-axis is the number of objects per 
r — W2 color bin. The black line shows the full source distribution. For red WISE sources 
with Wl — W2 > 0.8 (e.g., AGN candidates), two branches are apparent, one with blue 
colors (r — W2 ~ 4.5) and a secondary branch with r — W2 > 6. We propose the latter are 
type-2 AGN candidates. 
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Fig. 18. — Cumulative surface density of high-redshift ULIRG candidates selected by red 
optical-to-H^4 color. The single bar with the right side arrow shows sources without SDSS 
counterparts. 



- 35 - 




Fig. 19. — Optical spectra of the high-redshift ULIRG candidates selected by r — WA > 14. 
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Fig. 19. — Continue. 
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Fig. 20. — The colors and brightness as a function of redshift for the 8 r — WA > 14 ULIRG 
candidates with spectroscopic redshifts. The three panels are intended to investigate trends 
for selecting high-redshift ULIRG candidates. It is clear that fainter optical magnitudes or 
redder r — WA colors seem to select ULIRGs at higher redshifts. 
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